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ABSTRACT 


This report is an elimination of special nonlinearities of 
the Jeffcott equations in rotor dynamics . The immediate 
application of this analysis is directed toward understanding the 
excessive vibrations recorded in the LOX pump of the SSME during 
hot firing ground testing. 

Deadband, side force and rubbing are three possible sources 
of inducing nonlinearity in'the Jeffcott equations. The present 
analysis initially reduces these problems to the same 
mathematical description. A special frequency, named the 
nonlinear natural frequency, is defined and used to develop the 
solutions of the nonlinear Jeffcott equations as singular 
asymptotic expansions. This nonlinear natural frequency, which 
is the ratio of the cross-stiffness and the damping, plays a 
major role in determining response frequencies. 

Numerical solutions are included for comparison with the 
analysis. Also nonlinear f requenc y-reponse tables are made for a 
typical range of values. 

Finally the regions of stability of these nonlinear problems 
are bounded using inequalities of the system parameters. 
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1 . Introduction 

Vibrations are inherent in rotating machinery, Mathematical 
explanations of vibrations began with Jeffcott's description of 
the shaft's natural frequency of lateral vibration [5]. 
Unfortunately, Jeffcottfs linear model cannot account for all 
frequencies that have been observed experimentally. In 
particular, destructive vibrations have occurred in hot firing 
ground testing of the LOX pump of the SSME with no clue to these 
vibration's origins being offered by the linear model. 
Specifically, examination of the power spectral density (PSD) 
plots reveals unaccountable frequencies. Consequently, numerous 
investigations have been undertaken to study such rotors and to 
provide descriptions of the solutions of the two, coupled, 
second-order differential equations which describe the motion of 
the rotor's center of mass. Following the early work in 
rotordynamics by Yomarooto [7], one introduces a nonlinearity to 
the Jeffcott equations by including the effect of bearing 
clearance or deadband. In the pump, this deadband refers to the 
load carriers (ball bearings) and physically describes the 
clearance between the outer race of the bearing and the support 
housing. The work of Yomamoto did not include cross-stiffness, 
but a straight-forward derivation with this modification is 
easily obtained. A more limiting gap in his work is the 
assumption that the response is simply a perturbation of the 
forcing function. This is tantamount to assuming that one always 
has the graph of a circle as the solution. It is shown in this 
report that this generally is not the case, Both empirical 
results by Childs [1,2] and Gupta et al. [4] and numerical 
solutions using the fourth order Runge-Kutta algorithm by Control 
Dynamics Company [3] have been helpful in understanding the 
rotons motion for the nonlinear problem. This report extends 
the earlier work by using analytic expressions obtained from 
singular asymptotic expansions (method of multiple scales) to 
quantize the solution. 
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The primary objective of this report is to describe 
analytically solutions of the nonlinear Jeffcott equations with 
deadband (or an equivalent variation such as rubbing). To this 
end emphasis is placed on determining which frequencies one 
expects from the nonlinear problem, how these frequencies are 
related to the parameters of the differential equations, and 
where the solutions are stable. 

Section II is a heuristic development of responses to a 
nonlinear system. This discussion is based on a numerical 
solution and is intended to help develop some intuition about the 
solution of the nonlinear problem. 

Section III is the formal mathematical development of 
analytic solutions. This section begins with the discovery of 
the nonlinear natural frequency and proceeds to incorporate it in 
an asymptotic expansion of the solution. This leads to 
considering the two cases of resonant and non resonant driving 
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T frequencies. Herein also lies an explanation of why one expects 

1? the rotor's motion to be either a circle or an annulus, 

I * Section IV contains numerical examples which verify the 

} theoretical expansions. Typical frequency-response descriptions 
are also included in this pert. 

S Section V discusses bounds for stability regions. This is 

' done with comparison lemmas for majorizing the nonlinear Jeffcott 
equations . 

J Section VI concludes the report with a summary of the 

report's applications and directions for future studies. 
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2. Nondimensionalization and Heuristics 

The linear Jeffcott equations vhlch describe the displacement of the rotor 
center from its equilibrium position in the inertial, Cartesian coordinate 
system (y,z) (each in meters) are these: 


„ (1.) my-C g y-K g y-Q s z+mu(i) 2 coswc 

! (2.) oz«-(^ z+(^ y-K a z+rnuoi 2 sino) t 

where the shaft of the rotor lies along the x-axis and 


m ■ mass (kg.) 


, C * seal damping (kg./s.) 

* S 

„ K o seal stiffness (kg./s. 1 ) 

s 

1 Q « cross-coupling stiffness of seal (kg./s. 2 ) 

S f 

u a displacement of the shaft center of mass from the 
i geometric center (m.) 

« uj » angular velocity of the shaft (rad./s.) 

For the model to include bearing forces which hold the rotor in 
, position, one adds the terms • 

' ^/ )+u% //y 2 +z 2 ) 

and 

- yKg ( y-y6 / /f +z 2 )-K 0 ( z-z6//y *+ z 1 ), 
respectively, to right-hand sides of equations (1.) and (2.) 

f 

Here 

Kg » bearing stiffness (kg./s. 2 ) 

6 = clearance or deadband between housing and bearing (m.) 

y = coefficient of friction between housing and bearing (none). 

These bearing forces occur only when ✓y 2 +z 2 >6; otherwise, they are zero. 
Since y is nondimens lonal and typically small, one may regard it as zero 
without affecting the qualitative results. 

Equations (1.) - (2.) then become 

(3.) y + (C g /m)/+( 1/ m) ( +Kg ( 1-5/r ) ) y + ( Q s /m )z»u^ 2 cosw t 

(4.) z+(C s /m)z*-(Q s /m)y+( 1/m) [ K g +K 0 ( 1- 5/r ) ]z-ua3 2 sinwt 

when r » /j 1 +z 2 >6; otherwise, IL«0. Equations (3.) - (4.) can be 
put in nondlmensional form using a displacement g and a frequency 
o . One pair of candidates for g and o would be g «* 6 , the 

deadband size, and o 2 - u) 2 »K g +Kg the natural frequency of the 
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:t equations with deadband and mass 
lumerically using a fourth-order 
lese solutions are then plotted and 
lensity (PSD) investigation, the 
direction for initial, analytic 
isider the following special case: 

;-1.56 
1-200 ,000 . 
i- .0000285 

)=* 500Hz 

n after equations (1.) - (2.) and 
elative sense. Therefore, with 
00285 as the nond imensionali zing 
ons become 

/ Wj ) Z- 1 . 5 ( j/ujo ) 2 cos(u)/ U)o)T 
s 

(1-1/R) Z-1.5(w/w 0 ) 2 sin(w/ 

1/2 

y 2 +z 2 ) /6. Figure 1 shows the 

otion has become periodic for the 
wn . The graph of t vs. r = 6 R, 
at r has the form 








The prescribed initial conditions are used with the zero-order 
problem. The initial conditions for all higher-order problems ore 
all zer o . 

There is a circular argument involved here. One is using 
the answer, which was found numerically, to solve the problem. 
Although this is unacceptable mathematics, the results obtained 
in this heuristic discussion are not derived solutions but only 
suggestions of solutions. With this in mind, one proceeds to 
Bolve the zero - and first-order problems. The zero-order 
problem has solutions of the form: 

Y »C ( x) sin St +C (t >cos St +C cos<J>t -C sin$T 

fl 1 2 5 6 

Z »C (-r)sinST +C (t )cos St +C« sin<J>T +C, c.os 'J'T 

0 J w 3 

where C l( C 2 , Cj and C„ are all functions of time and all involve 
decaying exponential functions. Thus, in steady-state 

Y “C cos( + 4 > q ) 

3 ;ss i 0 

Z -C sin ( Af + A ) 

0 ,ss 7 v 0 

In the first-order problem, one is presented with the same 
left-hand side as the zero-order problem. Thus, in steady-state 
it is only the forcing function of the right-hand side that 
provides non-trivial solutions. Hence, one finds combinations of 
the form cos yx 'os <j>x and sin yt sin . 

Thus, the steady-state solutions of the first-order problem 
should look like these: 

Y l »C 9 cos( 4>i y ) T 

Z 1 «C 8 sin( 4>i y ) x 

Since higher order problems contain right-hand sides of the form 

(cos kyt)f(Y 0 , Z 0 , Y x , Z t Y k _ l , Z kl ) . 

it follows by induction on k that the kth-order steady-state 
solutions contain terms of the form 

Y c os ( ky ) x 

k 

Z ” sin( tp± ky ) t 

k 

Finally, one anticipates that the PSD of Y ( or equivalently Z) 
will show frequencies at <f>± ky . This is verified in Figure 3, 
which is the PSD of the (numerical) solution of the test problem. 
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3, General Theory 
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3*1 Nonlinear Natural Frequency 

Consider the homogeneous (E-0) equation corresponding to 
equation (7.). If this equotion were also linear (A*0), then 
exponentially growing or decaying solutions would generally 
result for a given set of system parameters. In the speciol case 
that (B/C) : -A+k , a sinusoidal solution is obtained with frequency 
B -B/C. To see this, consider the characteristic equation for 
W-exp(m t) : 

m 2 +Cm+{ A+k-iB ] «0 

m— C/2± {C 2 /4-A-k+i8 } 1/2 — C/2 ±{C 2 /4-( B/C ) 2 + iB } 1/2 

m— C/2* i {iC/2-B/C} 1/2 — C-iB/C, iB/C. 

In the nonlinear, homogeneous problem, k is replaced by k(l-A/R); 
hence, if R is a constant, then there is a wide spectrum of R for 
which (B/C) 2 may be A+k(l-A/R); i.e., if 

(10.) A<( B/C) 2 <A+k , 

then there is a constant value of R (with R>a ) for which 
(B/C) 2 »A+k( 1 -a/R) , This value of R is denoted by a and the 
corresponding frequency by 6 0 -B/C. This frequency is labeled 
the nonlinear natural frequency. Thus, whenever 

inequality (10.) is satisfied, equations (5.) - (6.) with E«0 
have steady-state solutions Y»a cos(3 0 t ) and Z«a sin (S 0 t). 

The same result 4 ' can be found from the polar equations (8.) 
- (9.) by assuming R and O' are constants. Then with E-0 those 
equations become 

(A+k(l-A/R)-(0') 2 ]R-0, Co’-B. 

Thus, 2 0 »G’“B/C and a-R-kA/ ( A+k- 3 fl 2 ) . 

Notice that g 0 -B/C<(A+k) =uo , the dimensionless natural 
frequency of the linear system. Thus, in considering the general 
nonhomogeneous problem, it is necessary to be aware of these three 
dimensionless frequencies: 

3 0 - the nonlinear natural frequency, 

.j 0 - the natural frequency, 

b - the driving frequency. 

Either 3 0 or u) 0 is an appropriate choice for o , the 
nond imensionalizing frequency. Correspondingly, one would select 
either a (with Bo ) or u (with u> B ) as the base displacement g. 
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3.3 Method of Multiple Scales. Op POOR QUALITY 

Nov consider a more formal development of analytic 
solutions. Straight-forward asymptotic expansions are not 
general enough for this problem since they always lead to a zero- 
order approximation of the form: 

W«M exp(iBoT ) + N exp(i 4>x ). 


The heuristic example considered in section 2 suggests that 
B 0 t must be replaced by 3 ( t ) which is close to but not exactly 
B 0 T • There are several methods of attacking such singular 

problems;, these .are described in Nuyfeh [61., Ope method, the 
method of averaging, is appropriate tor tnd JefFcott equations 

since it begins with the assupmtion that W = M (t ) exp (13 (x)) + 
N exp (1 (px ) . Another method, multiple scales, is also appropriate 
because one can envision the action of the rotor being based on 
two different time scales. The following development uses the 
method of multiple scales. The results are identical for the 
method of averaging. In this subsection the forcing frequency 
$ is restricted to differ from the nonlinear frequency g ; i.e., 
the problem is not resonant. Resonance, g ■ $ , is dealt with in 
the following subsection. 

Instead of one time scale x , assume the problem depends on 
many time scales: 


T 0 - x. T t - ex, T 2 «e z t , ... 

Henceforth, only T and T are used. Let W(t )*W(T o , T ) «W 0 (T „ , ) 

+-eWi (T , T 1 )+... Equation (11.) becomes a partial differential 

equation since 


d/dx -(3/3T )(d^ /dx ) + (3/9 T 1 )(dT l /d x)»D 0 +cD i 
and (d 2 /dt 2 ) = (D 4C D ) 2 . 

0 i 

Thus, one finds 

(12.) (D +£ D f (W + eW +. . . )+C(D + £D )(W ,+£W +...) 

0 1 o 1 u 1 J I 

+ < (W 0 +eW i + . . . )«ef (W fl +cW i + . . . ) + E<}> 2 exp( 14» T^ ) . 
Equating like powers of t yields 
(13.) e° :D 0 2 W 0 +CD 0 W 0 + kW 0 - E4) 2 exp(i<j)T 0 ). 

This is a linear problem with this steady-state solution 


W o =M exp( i6 0 T Q )+N exp(i<f>T 0 ) 

where N= E(p V ( - <J> 2 + i C <$)+< ) and M»M(T V ). To determine M one must 
examine the c -order problem and choose M to eliminate secular 
terms; see Nayfeh [6]: 
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c 1 •' D 0 2 V i+ CW i +<W i -2D 0 D l W 0 -CD ; W, + f(V 0 ). OH POOR QUAUTW 

With L-kA/t , the right-hand side of the last equation becomes 

-213 0 M’exp(iB 0 T 0 ) -CM ' exp ( i B 0 To ) 

+ L(1/|V 0 |-1) (M exp( i BoTo ) + N ex p( id> T 0 ) 1 

where I W 0 1-(1M| 2 + |N| 2 +MN exp[i(4>-3 0 )T 0 ]+MN exp[ i(0 o - )T 0 ] > 1/2 

To avoid secular terms one requires that the collective 
coefficient of exp(iB # T 0 ) be zero. Although an analytic solution 
of the differential equation f 0 r M(T X ) has not been found, one 
can qualitatively assess M based on a similar problem (van der 
Pol's equation) and specific numerical results (presented in the 
next section). 

Since M(T i ) is complex, it may be written as 

M(T 1 )«p(T l )exp[iS(T; )). Thus, 

A A 

w o“P(T i)exp[ iBo T 0+16 (T j) ]+N exp(i4>T 0 ) or, assuming S ( Tx ) is 
analytic near t«0, -p (T j) e x p [ i (Bo 4 cBi ) x + • • • 1 + N exp(i.^t). Thus 
the fundamental frequency of the nonlinear problem is not g 5 but. 

8 - B 0 + i -B l + . . . ; however, 0 must reduce to g Q when Ed> 2 °0* This 
frequency shift can account for the phenomenon of "tracking" that 
has been observed experimentally [2], Similarly, the frequency 
Y ° $ -S 0 that appears in. the expression for |W„| should be 
considered as Y =1-6. Then 1/|W 8 | shows all frequencies ny and 
W 0 /| W I shows all frequencies ny»£ , for n-0,1,... This suggests 
that M has a complex Fourier series of the form: 


n« -I s n exp(inyTi). 

Another factor of M must also be included since numerical 
examples show that M=0 if E<}> 2 is greater than some fixed value. 
This is sirilar to the behavior of the van der Pol oscillator; 
see [6]. Thus, one may speculate that M has a factor of the form 
F=l/ [ l+exp(- n T )] where n = n( EO 2 ) . This would imply that F ■* 1 
as T -* 00 when n >0 and F -* , 0 as T^^when n <0. Thus, M looks like: 


1 / [ 1 + ex p ( - nT r ) ] E s n exp{ inyT^ ) . 

n=-°° 

PSD plots of R show only frequencies of ny while plots of Y show 
frequencies of <5 and n yiB. 

3 . 4 Resonance . 

The resonance problem 3 can be handled in a manner quite 
similar to the nonresonance case. The approximation 4>=l+uE is 
made where E is the nond imens i ona 1 mass causing the imbalance and 
B 0 is the nond imensionalizing frequency. With the same multiple 
scale expansion as before, one finds the following two low-order 
problem: 
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with the solution W # »m(T^ )exp( IT 9 ) , OF POOR QUALTW 

0(E): s 2 W / aT 2 +C( 3 W / aT ) + [ 1- 1/a-iB )W + 23 2 W 0 / aT aT + 

C( 0 W o / IT )■( 1/E)( 1/ |W 0 j-l/a)W o + exp(iT # ) exp(i 0 T . 

Elimination of the secular terms of the latter equation requires 

0 — 2i M'exp(iT )-CM* exp( iT n ) +( 1/E )( 1 /M- 1/a )M exp(iT ) + 

e x p ( io T J exp(iT 0 ) 

or with M - mCTj) exp(i a (T t )) 

(2i+C)(m'+ici'ts)si(ni/E)(l/m-l/a)+ exp(ioTj) exp(-ia ). 

Now let ^ = o T t - Ci , separate real and imaginary parts, square 
both sides of the two resulting equations and add to find the 

steady-state solution (m'»0«a')s 

m 2 (l/a 2 -4 o /a + o 2 (c? +4 ) )+m(-2/a+4 a )-0. 

This magnitude result, a quadratic equation, is identical to that 
of Yomemoto [7] although the expansion there is with the base 
exp( i $ T 0 ) . 

Worn exp(i(T fl + a )) 

where m is found above and the frequency is 
(d/dt ) (T q + a ) = (d/dt )( t+ o Et- ip ) . 

But d i])/dt»0. Thus, the frequency is 1+ o E« q> . This phenomenon 
is often referred to as entrainment; i.e., the nonlinear 
frequency, G, is entrained by the forcing frequency, u) . 
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4. Examples . OF. POOR QUALTTK 

This section contains three numerical examples of the theory 
presented in Section 3. The first example shows how the solution 
changes for a fixed forcing frequency but increasing forcing 
magnitudes. Example 2 is similar to the first example but with a 
different set of parameters. Example 3 is a frequency-response 
illustration of variations in the forcing frequency and 

corresponding responses. 

Example 1. In this example the system constants used are 
these: y -0, m-1 lb.-s. 2 /in,, C_ = 240 lb.-s./in., K =0., 

K -1, 305,000 lb. /in., 0 -200,000 lb.7in., 6 -.0000285 in., and 

uj-500 Hertz-1000 it raa./s. Thus, g Q -833.33 rad./s. and 

a-. 000060913 in. The system is made nondimensional using a for 
the g-d isplacement and g for the o -frequency. With these 

choices, the constants of this equation: 

W"+CW* + [k(l- A / |W| )-iB]W=E0 2 exp(i «rr) 

have these values: C-.288, k-1.8792. A -.467865, B-.288, and 

<J>-6 1^5 . 

Figures 4 and 5 show changes in the solution Y vs. Z as E assumes 

the values 100n/( 1000 a for n-0,l,...,7. The graphs are 

plotted for ,2<t<.5s. The initial circle (for E-0) opens into Bn 
annular region, which becomes larger and thicker as E increases 
until a (transition) value of E occurs and the coefficient of 
exp(iSt) becomes zero. Thus, W-N exp(i4>t), a circle of radius 
|Nj. As E increased beyond this transition value, the solution 
remains a circle (figure 2.d) with radius |N| = | E<f> 2 / ( - + 10 $ 

+ k ( 1 — A/ 1 N l ) — i B ) | . 

The nonlinear natural frequency So is the angular 

frequency of the circle when E-0, but this frequency increases as 

E increases. This tracking phenomenon is displayed in the PSD 
plots of figures 6 and 7 using the dimensional frequency g . In 
these figures only the 120-180 Hertz range is shown. At 

E-7/ 10 , 000 tt 2 a , there is no frequency in this range; instead, the 
circle is tranversed at the forcing frequency, w . 

Figure 8, a typical full PSD plot, is the case 

E-4/10 , 000tt 2 a. As shown earlier, one expects frequencies of 

oj and g to appear, as well as harmonics of ny+g where y= -g 
and n=l, 2,... Thus, with g -150 Hertz, and oj«500 Hertz, one 
predicts that the PSD plot will exhibit peaks at 150, 200, 500, 

550, 850, 900, .. .Hertz. Figure 8 confirms these predictions. 

Example 2. In this example, the system constants from 

CDC[3] are these: y -0, m-. 20422 lb.-s. 2 /in., C =200 lb.-s./in., 

K =200,000 lb. /in., K. =1,000,000 lb. /in., u> =500 Hertz, 

Qg =C w/ 2 lb. /in., and 6-.0005 in. Thus, $=250 Hertz = 
500rr s rad . /s . and a». 0007183 in. Figure 9 summarizes changes in 
Y vs. Z as E varies by 250n/( 1000 it) 2 a, n=0, 1,...,4. These 
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t POOR QUALITY 

In general, if the time interval Is sufficiently large, the 
plots of annular regions will be completely filled (visually if 
not mathematically). However, when the ratio of 0 to y is p/q 
for small, positive integers p and q, then the curve actually 
falls on itself os time evolves and some attractive patterns 
appear. Figure 9. d is a case where 0/y=5/4, and the picture 
would look essentially the same whether shown for .45<t<.5s. 
(actual) or for .45<t<50s. On the other hand, figure 6.c is”more 
typical and would show a black annulus for ,45<t<50s. 

Example 3. 


As described above, one expects either a circle or an 
annulus at each fixed forcing frequency 4> • Furthermore, near 
resonance one always finds a circle since the characteristic 
frequency is entrained by the forcing frequency. Considering 
that the forcing function is a continuous function of this 
frequency $ , it follows that the solution is a circle or an 
annulus over an interval of the <f> -axis and that the number of 
transitions between the two shapes is finite over a bounded 
portion of the positive <(>-axis. 

TABLES I, II, and III illustrate these results for the 
problem with these parameter values: 


C “240 . / Wo 
s 

Q =200,000 . / w’. ; 
S ® 

K =0. 


5 =.0000285 
F =0.5, 1.0, or 1.5 
K =1,305,000. = w z 


Tables I, II, and III are for the three cases E=0.5, 1.0, 1.5, 

respectively. In all three cases, w 0 is used as the 

nondlmensionalizing frequency. Each of the three tables 

describes the response curve by listing the radius of circles or 
the inner and outer radii for annuli. Furthermore, the frequency 
y “ 1 G— uj | is given for annuli. The values listed for $ are 
every 0.1 except when 

1. more refinement is required to bound better the 

transition point 

2. less refinement is needed because of small variation in 
the response at different frequencies. 


Particular attention should be paid to Table III since this table 
shows three different regions in which the circle appears as the 
solution. The first two .tables contain only one such region. 
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5. Stability, 

To examine questions of stability for the nonlinear Jeffcott 
equations, it is helpful to use the following complex form: 

w+C g w+( Kg +K g ( 1- 6/r)-i(^ )w-F( u,u) 

where v*y+iz. Since the solution's stability is determined by 
the corresponding homogeneous equation, one has 

tf+C s w+[Kg+ K g (l- 6/r)-iQ g Jw-0. 

Furthermore, if a solution becomes unbounded then r > r fl for all 
constants r„ • Thus, as a comparison one may consider 

w+C g w+[K B +K g (U 6/r fl MQjw-O. 

Solutions of this equation are of the form w=exp(Mt). Hence: 

m 2 +C m+[K +K (1- fi/r )-iQ 3=0. 
s B s os 

The stability question is reduced to determining whether Re(m)<0. 
But 

m»-C / 2 ±[(C / 2 ) 2 — K -K (1- 6/r ) + iQ ] 1/2 
s s B s os 

and 

Re(m)=-C g /2 ±Re[C g 2 /4-Kg -K g ( 1-6 /r ) + i(^] 1/2 

*-C / 2 ± Re[C +iC ] 1/2 
s 12 

where 

C =C 2 /4-K -K (1- 5 /r ) 
is B s a 

VI • 

In this modified form, the stability question becomes this: 

1/2 

Solutions are stable if and only if ReCC^iQ ) <C S / 2 . 

The following lemmas are useful for bounding the value of Q s /C s 
in order to find stable solutions. 

Lemma 1. Re( C 1+ iC 0 1 /2 =(|C X I / 2) 1/2 [ ( l + ( C /C ) 2 J 1 /2 ± 1 1 1 /2 vher e the 
minus sign is used for C i < 0. 

Proof: Case C < 0. C < 0 implies j C, ! =-C . Then: 

i i 1 1 * i 

Re(C +i.C ) 1/2 
i 2 

»(C 2 +C 2 ) 1 / 4 cos[l/2 Atan (C, /- |C , | ) ] 
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tj. 


«(C 2 +C 2 ) 1/4 c o s [ ( it — e)/2] 

1 2 

where 9 » Atan ( C ^ / 1 C \ ) 

-(Cj +C 2 ) 1/4 sin(0/2) 

-(C 2 +C 2 ) 1/4 [ (1-cos 0) / 2 ] 1/2 

-(CJ +C 2 ) 1/4 (1/2) 1/2 [1-ICJ/CCJ +C 2 ) 1/2 J 1/2 

-1/2 1/2 [(C J +C 2 2 ) 1/2 - |C 1 | ] 1/2 

»(|cj / 2 ) 1/2 [(1 + (C 2 /C x ) 2 ) 1/2 -1J 1/2 


Case c >0. C >0 implies C *|C | . Then 


l l 


Re(C +iC J /2 
1 2 

“(C 2 +C 2 ) 1A cos[1/2 AtanC^/jCj)] 

= (C 2 +C 2 )1A[(1 + cos Atan C /I C |)/2]1/2 
12 2 1 

-cc* +c 2 ) 1/4 [ C 1+ |c i I /(C 2 +C 2 ) 1/2 )/2] 1/2 

12 1 1 2 


*1/2 1/2 [(C 2 +C 2 ) 1/2 + |CJ ] 1/2 
*( |CJ /2) 1/2 [ ( 1 + (C 2 /C i ) 2 )+l] 1/2 


Lemma 


2- If (Q s /C s ) 2 <1^ f then for all & _<r o < °° t 


Re(C +iC ,) 1/2 <C /2 . 

12 g 


Proof: Case C ^<0 . If C y <0 , then |C j «-C • 

(Q /C ) 2 <K 


-K <-(Q /C ) 2 
B s s 



-(K +K (1- 6 / r )) < -K < -(Q /C ) 2 
SB o B s s 

C * ( C / 2 ) 2 -(K +K(1- 6 /r ) ) < (C /2) 2 -(Q /C ) 2 

IS SB 0 S SS 

-ic,i < 

-*V PJ < c; -4Qj 

4Q 2 < C ‘ +4C 2 |C I 
s s s 1 1 1 

4(Q 2 +C 2 ) < C “ +4C 2 |C | +4 |C ! 2 

s 1 S S 1 1 1 1 1 1 

4(Q 2+C, 2 ) < (C 2 +2| C x | ) 2 

5 S 

2|Cx| [U(Q /C,) 2 ] 1/2 < C 2 +2|Cx| 
s s 

{1 + (Q /Cj ) 2 ] 1/2 < C 2 /2 1 Cl j +1 


|C |/ 2 <[l+<Q a /Ci)M 


1/2 


-1) < C 2 /4 
s 
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( lc x I /2) 1/2 ( [ 1 + CQ 0 /C r ) 2 ] 1/2 -U < c /2. 


But by Lemma 1, the last left-hand side is Re(Ci+iCj )^ 2 . Thus, 

Re(C. +iC s ft 2 < C 12. 

11 a 

Case C 1 >0, This case is identical to the previous argument 
except that one has + | C x | instead of -|Cj on the fifth line of 
the inequalities and thereafter the corresponding term switches 
signs . 


Lemma 3. If ( Q / C ) 2 > K n +K , then for all <5 < r <® 

SS D S — » 0 

ReCCi + iCj)^ 2 > C J2. 

Proof: After noting 


<V C S > ! > V K S > %+V»- « /r 0 > 

one has 


-k b -k a- { / r t ) > -(lyc^ 


and 


C x — K -K (1- 6/r 0 ) + ( C 2 / 2 ) 2 > (C Jl ) 2 -(Q JC ) 2 

Do S 5 5 

The remaining steps in the derivation are identical to Lemma 2 
for both cases of, C, except that the inequalities are reversed. 
Thus, Re(C 1 +1C z y ^ > C g /2. 

Figure 10 summarizes the stable and unstable regions as 
functions of (Q /C ) 2 for the linear and nonlinear homogeneous 
problems. In this figure, there is only one point, on the 

(Q„/C g ) 2 -axis at which oscillations may occur for the linear 
problem. For the nonlinear case, oscillation and hence 

transition between stable and unstable still occurs at only one 
point on the (Cl / C g ) 2 -axis , but that point may be anywhere 
between Kg and +k , 
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transition point to consider. These problems, however, greatly 
extend the model’s mimcry of an observed rotor's behavior. 

This report is also limited to single forcing functions. 
Realistically, one must consider multiple forcing functions. 
Here again, superposition will fail for the nonlinear problem 
although it may be a first approximation. Certainly, harmonics 
(oums and differences of driving frequencies) may appear end 
could possible dominate. These problems appear to introduce no 
new theory, but do increase the computational complications. 

Stability for all these problems remains the central focus. 
Even in the symmetric nonlinear problem with a single driver, it 
is still an open question of whether the response may move from 
an annulus to a circle (or vice versa) when it is perturbed. 




7 , References : 

1. Childs, D. W., '’The Space Shuttle Main Engine High- 

Pressure Fuel Turbopump Rotordynamics Instability 
Problem", Trans. ASME, Journal of Engineering for 
Power . Jan. 1978, pp. 48-57. 

2. Childs, D. W., "Rotordy namic Characteristics of the 

HPOTP (High Pressure Oxygen Turbopump) of the SSME 
(Space Shuttle Main Engine)”, NASA MSFC Contract 
NAS8-34505, Turboraachinery Laboratories Report 
RD-1-84, 30 January 1984, 

3. Control Dynamics Company, "Effects of Bearing Deadbands 

on Bearing Loads and Rotor Instability”, NASA MSFC 
Contract NAS8-35050, 20 January 1984. 

4. Gupta, P. K., Winn, L. W, , and Wilcock, D. F., 

"Vibrational Characteristics of Ball Bearings", 
Journal of Lubrication Technology, ASME Trans., 

tov^f ,“No:“"r9rr;“pp‘.‘ 

5. Jeffcott, H. H. , "The Lateral Vibration of Loaded Shafts 

in the Neighborhood of a Whirling Speed - The Effect 
of Want of Balance", Philosophical Ma gazine , Series 
6, Vol . 37, 1919, p. 3 W. 

6. Nayfeh, A. H., Perturbation Methods, J. Wiley & Sons, 

1973. 

7. Yamamoto, T. T., "On Critical Speeds of a Shaft”, 

Memories of the Facul t y of Engineering , Nagoya 
(Japan) University, Vol.~ 5, No” IT54, 


19 









X10 


Mm 




ORIGINAL PAGE to 
OE POOR QUALITY 







=} o 
o 

U_ E 

i_ 

4— O 
O 4— 

00 C 
- C — • 
CTO'-' 











V V. 


r* 

\ H 

k 




(a.) 


ORIGINAL PAGE IS 
Of- POOR QUALITY. 









ORIGINAL PAGE & 
( b .) OR POOR QUALITY 













ORIGINAL PAGE IS 
OF. POOR QUALITY 


Figure 9 
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ORIGINAL PAGE tS 
TABLE I ° E P0 °R QUALITY! 

E - 0.5 


4 > 

Shape 

Radius 

(radii) 

y (Hz . ) 

(Circle/Annulus) 




0.0 

C 

2.137 

— 

0.1 

A 

2.130 - 

2.145 

115 

0.2 

A 

2.105 - 

2.169 

97 

0.4 

A 

1.959 - 

2.305 

61 

0.5 

A 

1.747 - 

2.470 

41 

0.53 

A 

1.618 - 

2.549 

34 

0.54 

C 

1.597 

_ 

0.6 

c 

1.833 

- 

0.7 

c 

2.441 

- 

0.8 

c 

3.656 

- 

0.9 

c 

6.954 

- 

0.92 

c 

8.249 

- 

0.93 

A 

1.212 - 

4.018 

13 

1.0 

A 

1.369 - 

3.230 

40 

1.2 

A 

1.551 - 

2.858 

81 - 

1.4 

A 

1.643 - 

2.822 

118 *' 

1.6 

A 

1.648 - 

2,748 

154 

1 .8 

A 

1.694 - 

2.633 

190 

2.0 

A 

1.804 - 

2.679 

222 

2.2 

A 

1.939 - 

2.830 

263 

2.4 

A 

1.993 - 

3.162 

294 

2.6 

A 

1.632 - 

3.125 

333 

2.8 

A 

1.533 - 

2.864 

370 

3.0 

A 

1.528 - 

2.769 

410 


r 


G 


30 
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F. - 1.0 

Shape Radius (radii) y( 

(Circle/ Annulus) 


0.0 C 

0.2 A 

0.4 A 

0.41 A 

0.42 C 

0.5 C 

0.6 C 

0.8 C 

0.9 C 

0.97 C 

0.99 C 

1.0 A 

1.2 A 

1.4 A 

1.6 A 

1.8 A 

2.0 A 

2.2 A 


2.137 

2.073 - 2.201 
1.767 - 2.463 
1.737 - 2.486 
1.398 
1.650 
2.118 
4.548 
9.174 
18.524 
10.436 
.87: 

1 .oo: 

1.14, 

1.43 
1.15 
1.34 
1 . 591 
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1 a b le 1 1 1 0F pQ0R Q UAL mr. 

E ■ 1.5 


* 

Shape 

(Circle/ Annulus) 

Radius (radii) 

y (Hz 

0.0 

C 

2.137 


0.2 

A 

2.040 - 2.230 

96 

0,3 

A 

1.885 - 2.371 

78 

0.35 

A 

1.755 - 2.476 

68 

0.36 

C 

1.364 

* 

0 . A 

C 

1.449 

- 

0.6 

C 

2.398 

* 

0.8 

C 

5.427 

- 

0.9 

C 

11.269 

- 

0.98 

C 

27.247 

- 

0.99 

c 

25.102 


1.0 

c 

19.448 

- 

1.2 

c 

2.585 

— 

1.38 

c 

2.028 

* 

1.39 

A 

.357 - 4.628 

99 

1.4 

A 

.529 - 4.832 

100 

1.6 

A 

.837 - 4.660 

139 

1.77 

A 

.359 - 3.832 

172 

1.78 

C 

1 .717 


1.8 

C 

1.709 

— 

1.91 

C 

1.676 

— 

1.92 

A 

.100 - 3.121 

195 

2.0 

A 

.496 - 3.282 

217 

2.2 

A 

1.062 - 4.076 

250 

2.4 

A 

1.496 - 4.564 

286 

2.6 

A 

1.692 - 5.341 

312 

2.8 

A 

.222 - 3.987 

359 

2.81 

C 

1.565 

— 

2.9 

C 

1.560 


3.0 

C 

1.555 

— 


